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Abstract 



> 

O ' The dynamics of passive Brownian tracer particles in steady two-dimensional potential 



flows between sources and sinks is investigated. The first-passage probability, pit), 



^ I exhibits power-law decay with a velocity-dependent exponent in radial flow and an 

I order-dependent exponent in multipolar flows. For the latter, there also occur diffusive 

^ ' "echo" shoulders and exponential decays associated with stagnation points in the flow. 

I For spatially extended dipole sinks, the spatial distribution of the collected tracer is 

S ' independent of the overall magnitude of the flow field. 

^ ■ P.A.C.S. Numbers: 47.55.Mh, 05.40. +j 

o ■ 
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The motion of passive tracer in a flowing fluid is a convenient diagnostic tool for moni- 



H ' toring the characteristics of the flow field, as well as properties of the background medium. 

In porous media flows, particularly in groundwater or hydrocarbon recovery, tracer measure- 
ments often are the only source of internal information about the system [1]. Considerable 
effort has focused on one- dimensional flows where tracer is introduced at one end of the 
sample and the transit time distribution is monitored upon exit at the other end. However, 
in many situations, such as fluid pumped into an "injection" well and extracted from one or 
more "producing" wells, the flow is more likely to be radial, or multipolar, rather than linear. 
Our goal, in this letter, is to provide general insights about the first-passage properties for 
dynamically- neutral tracer in such flows [2]. 

We shall focus on steady two-dimensional flows which are generated by multipolar con- 
figurations of sources and sinks. These appear to encompass many cases of possible physical 
relevance, as well as illustrating a rich range of phenomenology. Mathematically, the first 



passage of passive Brownian tracer particles is determined by their concentration, c(r, i), 
which satisfies the convection-diffusion equation (CDE), 

+ u{r) ■ Vc(r, t) = DV'c{f, t), (1) 

with the velocity field u{r) arising from a two-dimensional potential. We will elucidate 
the features of p{t), the distribution of transit times or the first-passage probability between 
source and sink, to characterize the motion of dynamically-neutral tracer in steady fiows. 

First consider radial or monopolar flow, u{r) = Uor/r. Non- universal features arise due 
to the competition between the flow and the centrifugal term in the Laplacian, which both 
vary as 1/r. Consequently, the form of the CDE in two dimensions {d = 2) with uq finite is 
the same as that of a system with no drift but with d ^ 2. Thus changing the amplitude 
uo is equivalent to changing the spatial dimension, leading to non-universal first-passage 
properties as a function of the drift. 

To illustrate these features, we have solved for p{t) in a circular geometry with an inner 
absorbing radius a, infinite outer radius, and an initial ring of tracer particles which are 
released at r = Tq. This first passage probability coincides with the radial fiux at r = a, i.e., 
p{t) — —Ddc{r,t)/dr + u(r)c(r,t). Using standard Green's function methods, the Laplace 
transform, p{s) = p{t) e~^^ dt, is [2] 



.^yX^^o^^ (2) 

where Ky is the modified Bessel function of the second kind, and v = Uq/2D = Pe /2, 
where Pe is the Peclet number. The long-time behavior is determined from the small-s 

expansion of p, which generically has the form, p(s) ~ p(0) — as"" -\-bs^ -\- By construction, 

the leading term p{s = 0) = p{t) dt = E{ro) is just the probability that a tracer particle 
which starts at r = ro eventually reaches the absorbing circle. Time-dependent properties can 
be inferred from the correction terms in Eq. (2). When a < 1, then the mean first-passage 
time to the absorbing circle, (t), diverges and the first-passage probability in real time has a 
power law tail which varies as t~^"~^^\ However, when a reaches unity, the coefficient of —s 
equals (t). 

For outward fiow, the small-s expansion of p{s) is 

Thus E{ro) = (a/ro)'^'' < 1 and some tracer escapes to infinity. This reflects the equivalence 
to a purely diffusive system of spatial dimension d > 2, for which the eventual return 
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probability is less than unity [3]. Prom the correction term, we deduce that (t) — > oo for 
< u < 1, while for i/ > 1, {t) is finite. Since E{ro) < 1, {t) in the case u > is conditional 
on those particles which actually hit the absorber. Note that for sufficiently large outflow, a 
tracer particle must reach the absorber in a finite time if it is to be absorbed at all. 

For inflow {u < 0), the first passage probability has the small-s expansion 

where /j, = —u > 0. Thus tracer hits the absorbing circle with certainty, while p{t) is non- 
universal and varies as t~^^'^^\ Additionally for fi < 1, {t) is infinite, while for // > 1, {t) is 
finite. Even though the tracer is guaranteed to reach the absorber for any inflow (including 
zero), the first-passage time is finite only if the drift is sufficiently strong. 

We now discuss tracer motion in a fully two-dimensional dipole flow field due to a point 
source and sink at (0, =Fa) within a circle of radius R, with total flux ±Q. The velocity field 
may be written as 

^ ^, ,dilj dij 

Here the velocity potential (j) and streamfunction ip are the real and imaginary parts of the 
complex potential 

^ = + ^V; = ^ln) 6 

ZTT [z — a){z — R^/a) 

where ±i?^/ a are the image locations to impose vanishing normal velocity on the boundary, 
and z — X + iy. We assume that fluid and tracer are emitted uniformly in angle about the 
source. 

In the limit of pure convection in an infinite plane (R — > oo), the tracer motion is 
determined by the velocity field, df/dt — u{f{t)). Since the second equality in Eq. (5) 
represents a Hamiltonian system, we have applied a canonical transformation technique to 
compute particle trajectories analytically [4], and thereby determine that the transit time 
from source to sink for a particle emitted at polar angle 6 with respect to the source is 
t{6) = 2csc'^ 9 [1 — ^^cot^^ ]; see also [5]. Although the motion is deterministic, a distribution 
in transit times results from the distribution in initial angle. Since all streamlines lead from 
source to sink, there is a one-to-one equivalence between the angular and time dependence 
of the first-passage probability. This leads to 

S{2T:)-^/'-r'/\ (7) 



p{t)^p{e) 



dt 



a 3 ' 



due to tracer particles which initially move away from the sink on distant dipole streamlines. 
It is this aspect of the distribution which explores the global structure of the flow field. Note 
that the mean transit time, (t) = J^^ dttp{t), diverges. 
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To incorporate molecular diffusion, we resort to numerical simulation, and have employed 
several complementary methods. The simplest is grid-free Monte Carlo time-stepping us- 
ing individual random walkers, which works best at high Peclet numbers. In time At, a 
walker is displaced by Ar = u(r)At + n(4At/Pe)^/^, where n is a unit vector of random 
orientation. Figure 1(a) gives a typical result for Pe = 2700 with N = 500,000 random 
walkers. There is an early-time peak, followed by the anticipated t^^/^ power-law decay, and 
then a noisy exponential region dominated by diffusion. To clarify the latter domain, we 
have devised a lattice "probability propagation" algorithm. A probability element at 
translates through a distance Af— ■u(r)At in a time At to an off-lattice position {x, y), and 
is then redistributed among the 5-site nearest lattice neighborhood of the target position. 
The redistribution rule is chosen to ensure that the average displacement of the 5-site group 
remains equal to Ar and that diffusion is spatially constant by imposing a constant fluctua- 
tion in the displacement of this group. While this method has no statistical fluctuation, the 
simple form of the redistribution rule restricts the method to low flow rates. 

The probability propagation results are equivalent to standard finite-difference meth- 
ods, but programming is extremely simple. A typical result for p{t) (Fig. 1(b)), exhibits 
the following four generic features: an early-arrival regime, a power-law decay, a "diffusive 
echo" shoulder (whose resolution in simulations is dependent on the absence of statistical 
fluctuations), and an ultimate exponential decay. The shoulder stems from particles which 
reflect from the boundary before reaching the sink. Because of the low Peclet number in- 
volved, non-trivial asymptotic exponent estimation and extrapolation [2] are required to 
verify the exponent in the power-law region. Since p{t) ~ l/{t\n^t) in the diffusive limit [3] 
and p{t) ~ as Pe — > oo, it seems unlikely that a simulation method which cannot be 
applied at large Pe will be able to resolve these two hmiting behaviors cleanly. 

The decay at the longest times may be correlated with stagnation points in the flow 
field which generically lead to an exponential decay in the transit-time distribution [6]. 
Without loss of generality, consider a stagnation point at the origin with the local velocity 
u — df/dt — {—Gx,Gy). If the stagnation point is approached along the x-axis, then 
trajectories which pass near the origin and then escape have the form r(t) = [x^e'^^ ^y^e^^) 
with I/O -C Xq. The time T spent near the stagnation point can be defined by requiring 
that y ~ ?7o, where Uq is some characteristic 0(1) velocity. Since y ~ y^Ge^^ , this gives 
T ~ In (f/o/yoC). Typically, for trajectories which pass near a stagnation point, the 
transit time between source and sink is dominated by this value of T. The corresponding 
distribution of transit times may now be obtained by accounting for the distribution of initial 
positions yo- Consequently, 



p{T) ^Po{yo) 



dyo 



Po(ge-«^) 



-Ue 



-GT 



Po(0)e-«^. (8) 



dT 

Thus the transit time distribution depends on the local shear rate near the stagnation point. 
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but not on details of the initial spatial distribution of tracer po- In the present example, 
applying Eqs. (5), (6), and (8) to the stagnation points at (±i?, 0) gives an exponential decay 
rate with G oc R~^. The decay rate observed in simulations will typically be the smaller of 
this value and the diffusive decay rate, which is proportional to for the dipole case. 

For a source-sink of arbitrary multipolarity, the exponent of the power-law tail in p{t) 
depends on the multipole moment. This exponent can be obtained by the following simple 
argument. For a 2^-pole, the leading behavior of the complex potential is ^ ~ so that 
the corresponding streamlines arc il) = Im\& ~ r"''^ sin iV^^, for suitable orientation of the 
axes. In the convective limit, particle trajectories are defined hj ip = const.. This constraint 
implies that ilj{r,6) = ilj{R,'ir/2N) = R^^ , or = R'^ sin N6, where R is the maximum 
distance from the origin on the streamline. The angular velocity on such a streamline is 



^ ^ ^ Nr^'-'smNO ~ i?-^-^ sin'^/^ {N9). (9) 

dt rdr ^ ^ ^ ^ 

The transit time on this trajectory is the time required for 9 to vary between and ir/N. 

Using the above approximation for the angular velocity, the transit time scales as t — 

Jq^^ d6 (dO/dt)"^ ~ To determine the transit time distribution, we relate the value 

of R to the initial angle of emission 9q at the source. Suppose the source is at the origin 

and oriented so that the long-excursion streamlines arc associated with 9q 0. For r ^ e, 
= R~^ sin 9o ^~^9q, or 9o ~ R~^ . Thus the transit time probability distribution 

is 



p{t)=p{9o) 



d9o 


1 


d9Q dR 


dt 


" 2n 


'dR~dt 



While these predictions have been verified numerically in the convective limit, there are as- 
yet unexplained diffusive features in probability propagation simulations. For example, the 
form of pit) for a quadrupolc consisting of charges {—Q, 2Q, —Q) is quite different from that 
of an "inverted" quadrupole {Q,—2Q,Q). 

Another important diagnostic is the distribution of where on the sink the tracer is col- 
lected. In the dipole case, with passive tracer released in steady potential flow between a 
source and a single spatially extended equipotential sink, we demonstrate that the time- 
integrated tracer flux distribution at a given point on the sink is proportional to the in- 
coming fluid velocity at this point. By linearity, this implies that the spatial distribution of 
the collected tracer is independent of the overall magnitude of the flow. This establishes a 
useful general equivalence between tracer distributions for pure diffusion and in combined 
convection and diffusion. 

To prove the theorem, we use the complex potential as the conformal mapping to trans- 
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form the flow domain from {x, y) into (0, ■?/')■ In these variables the CDE becomes 

This transformation to linear flow introduces a spatially dependent time step which is sin- 
gular as (0, ivr) is approached. In the (0, ip) coordinate system, the source has co-ordinate 
= 00- Foi' the initial condition of a delta-function pulse of tracer injected at the source, then 
in the equation of motion for ^{4>,ijj) = Jl^ dtc{t,(l),ijj), the time derivative term integrates 
to zero, yielding 

with boundary conditions (1 — Pe~^9/90)^ = a Pe -independent constant, at = 0o, and 
^ = at = 01. The solution of Eq. (12) is ^ = K{1 - e-P^('^-<^i))^ and the time-integrated 
tracer flux arriving at angular position '0 on the sink, ^(-0), is 



P(^) ^ -Pe - /; ,t ^(0., ^, t) = -Pe (01, ^) = -Pe - 



(13) 

(01 .V") 



where h is the unit normal to the sink. In the last expression, the first factor is a constant, 
and the second is simply the normal velocity of fluid at the sink. Hence p('0) oc Un{ip) as 
claimed. 

Notice that the theorem also holds in the pure diffusion limit, as can be seen by taking 
Pe ^ either in the equation for C, or in its solution. In this case, the appropriate statement 
is that ^("0) is constant for a constant arclcngth of the sink. Conversely in the limit of 
no diffusion the theorem is obvious, because then tracer particles remain on their initial 
streamline, and the tracer flux is simply proportional to the fluid flux. The nontrivial 
implication is that the local integrated flux to a simply-connected sink is independent of the 
Peclet number, although the transient flux is Peclet number dependent. 

In summary, we have investigated the first-passage probability p{t) of dynamically-neutral 
tracer in two-dimensional potential fiows. For radial flow, convection and diffusion arc of 
the same order leading to non-universal flrst-passage characteristics. In general multipole 
flows, the exponent of p{t) depends on the multipole order. The effects of diffusion appear 
to be subdominant for the dipole but relevant for the quadrupole, as qualitative features of 
p{t) depend on the sense of the quadrupolar ffow when diffusion is present. The influence 
of stagnation points on the asymptotic properties of p{t) was determined. Finally, the time- 
integrated flux to a flxed arclength of sink is independent of the overall magnitude of the 
flow fleld. 
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Figure Caption 

Figure 1. First-passage probability distribution in combined dipolar flow and molecular 
diffusion within a circle of radius R = 400 with source and sink at a; = =f20, 
obtained by (a) following the motion of individual random walks for Pe — 2700 
and (b) probability propagation for Peclet numbers Pe — 1.6, 0.8, 0.4, and 0.2 
(tallest to shortest peak, respectively). The inset in (a) shows the streamhnes. 
Also in (a), the abscissa is the dimensionless time t ^ t ■ (27ra^/(5). 
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